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The experiment E97-006 was performed at Jefferson Lab to measure the momentum and energy distribution of 
protons in the nucleus far from the region of the (approximate) validity of the mean field description, i.e. at high 
momentum and energies. The occurrence of this strength is long known from occupation numbers less than one. 
In the experiment reported here this strength was directly measured for the first time. The results are compared 
to modern many-body theories. Further the transparency factor of 12 C was determined in the Q 2 -region of 0.6 
to 1.8 (GeV/c) 2 . 



1. INTRODUCTION 

The energy and momentum distribution of nu- 
cleons bound in the nucleus are obviously deter- 
mined by the nucleon-nucleon potential. The N- 
N potential consists of a repulsive core at small 
nucleon-nucleon distances of less than rs 1 fm and 
an attractive part. Such a realistic potential is 
usually constructed from N-N scattering data. In 
a vector meson exchange model the long-range 
part is represented by an exchange of a virtual 
pion whereas the short-range part is provided by 
heavier mesons like the omega and the rho. The 
strongly repulsive part, which is caused by the 
short-range and tensor components of the inter- 
action, is ignored in the mean field description. In 
the Independent Particle Shell Model (IPSM) it 
is assumed that the nucleons move independently 
from each other. The repulsive part responsible 
for N-N correlations is neglected. Therefore all 
nucleons reside below the Fermi energy and Fermi 
momentum of about 250 MeV/c. All orbits are 
fully occupied in the IPSM. It was somewhat of a 
surprise when (e,e'p) measurements performed at 
NIKHEF 1 found values of typically 65 % for the 
spectroscopic factor of the valence orbits in nuclei 
from lithium to lead, ie. 35 % are missing. Short- 
range (SRC) and tensor correlations are responsi- 
ble for up to 15 % depletion. Due to the repulsive 
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core of the interaction at small distances nucle- 
ons are scattered to high energy and momentum. 
Another 20 % is coming from long-range correla- 
tions (LRC) which lead to a fragmentation of the 
single-particle states due to coupling to collective 
modes like vibration. This influences the spectral 
function in particular at small energy, i.e. the 
valence states. For the deeper lying orbits the 
occupation is larger. 

Modern many-body theories like the Correlated 
Basis Function theory (CBF) [2], the Green's 
function approach (GF) |3J and the self-consistent 
Green's function theory (SCGF) @] are able to 
deal with the realistic N-N potential directly and 
therefore contain short-range and tensor correla- 
tions. The remarkable difference to the momen- 
tum and energy distribution derived in the IPSM 
is the additional strength found at large energy E 
and momentum k. Above the Fermi momentum 
Icf the single-particle strength is negligible. As 
a consequence the search for SRC has to concen- 
trate on large E and k. A practical tool is the 
(e,e'p) experiment. 

2. EXPERIMENT 

The experiment was performed in Hall C at 
Jefferson Lab employing three quasi-parallel and 
two perpendicular kinematics at a q > 1 (GeV/c) 
(for a detailed discussion see [5]). Electrons of 
3.3 GeV energy and beam currents up to 60 /iA 
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Figure 1. Coverage of the E m ,p m -plane by the 
runs taken in parallel kinematics shown in a cross 
section times phase space plot. (Due to the large 
momentum acceptance of the spectrometers, part 
of the data (green) are for 8k q > 45°). 



were incident upon 12 C, 27 Al, 56 Fe and 197 Au 
targets. The scattered electrons were detected 
in the HMS spectrometer (central momenta 2 - 
2.8 GeV/c), the protons were detected in the SOS 
spectrometer (central momenta 0.8 - 1.7 GeV/c). 
Fig. 2] gives the kinematical coverage for the par- 
allel kinematics. 

From the measured scattering angles and mo- 
menta of the electron and the proton the missing 
energy 



E m — E e — E e i —Tp> —Ta-i 
and the missing momentum 
Vrn = q-Pp' 



(1) 



(2) 



can be constructed. Here E e is the electron beam 
energy, E e > the energy of the scattered electron 
and Tpi (Ta-i) the kinetic energy of the knocked- 
out proton (the residual nucleus). In Plane Wave 
Impuls Approximation (PWIA) there is a direct 
relation between the measured quantities and the 
theory In this case E m can be identified with 
the removal energy E of the proton in the nucleus 



and p m with its initial momentum —k. However, 
in the nucleus the hit proton can interact with 
the other nucleons of the nucleus. Obviously this 
rescattering contribution is reduced for large mo- 
mentum transfer q and increases for heavier nu- 
clei. Calculations of the rescattering contribution 
[H] confirm that rescattering is small for 12 C in 
parallel kinematics where the momentum trans- 
fer q is parallel to the initial momentum of the 
proton. The data taken in perpendicular kine- 
matics lead to a three times larger strength in- 
tegrated over the experimental acceptance and 
compared to the parallel kinematics. Rescatter- 
ing contributions cannot entirely explain this dif- 
ference. Probably multi-step processes and Me- 
son Exchange Currents (MEC) need to be con- 
sidered. Work in this direction is underway |7I9| . 

The raw data were analyzed using two differ- 
ent procedures, both based on an iterative ap- 
proach and a model spectral function. In one, 
the phase space is taken from a Monte Carlo sim- 
ulation of the experiment, and the spectral func- 
tion is determined from the acceptance corrected 
cross sections. Radiative corrections are taken 
into account according to [S] ■ The other is based 
on a bin-by-bin comparison of experimental and 
Monte Carlo yield, where the Monte Carlo simu- 
lates the known radiative processes, multiple scat- 
tering and energy loss of the particles using spec- 
trometer transfer matrices. The parameters of 
the model spectral function then are iterated to 
get agreement between data and simulation. We 
have found good agreement between the two pro- 
cedures. 

In PWIA the spectral function can be obtained 
via the absolute differential cross section 



da 



dQ e dil p dE e > dEpi 



= K a ep S(E m ,pZi) Ta(Q 2 ).(3) 



Here K is a kinematical factor and a ep the e-p 
cross section for a moving proton bound in the 
nucleus. Since the bound proton is off-shell, i.e. 
the usual relation between mass, energy and mo- 
mentum is not longer valid, the expression for the 
e-p cross section cannot be derived uniquely. Dif- 
ferent versions are on the market like a cc \ and 
cr CC 2 10 which are derived from two expressions 
of the electromagnetic current operator T\ and 
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Figure 2. Nuclear transparency Ta for C, Fe and 
Au as a function of the proton kinetic energy T p 
compared to the correlated Glauber calculations 
(solid lines). The data indicated by circles are 
from the NE18-experiment at SLAC ^2jj squares 
and diamonds are Jlab data of ^Sj and an( i 
from Bates (triangle down). The result in- 
dicated by stars is obtained with the correlated 
spectral function of |16| . 



T2- Both are equivalent in the on-shcll case but 
differ for an off-shell proton. In the following the 
e-p cross section a cc |S] will be used where the 
kinematical variables are not modified to obtain 
an equivalent on-shell kinematics as in the case of 
a cc i and a CC 2- The difference between the cross 
sections is small in the IPSM region but increases 
with increasing p m and E m . The choice of a cc is 
motivated by the fact that it gives the best agree- 
ment of the spectral function for same E m and p m 
but measured in different (parallel) kinematics. 

The absorption of the proton in the nucleus 
is accounted for by the nuclear transparency Ta- 
Also for a check the analysis procedure Ta was 
measured for five different kinematics on 12 C. The 
results are presented in the next section and com- 
pared to the correlated Glauber theory of Benhar 

m. 

3. TRANSPARENCY 

Measurements of the transparency rely on com- 
parison of the yield measured for a specific nu- 



cleus to the yield calculated assuming PWIA. 
This transparency accounts for absorption of 
the proton in the nucleus but also for large 
momentum changes due to an interaction with 
the surrounding nucleons. The calculated yield 
N mm (E m ,p m ) is taken from a Monte Carlo sim- 
ulation which contains the detector acceptances, 
radiative corrections etc. It is checked that the 
spectra from the simulation are in good agree- 
ment with the one measured or reconstructed 
from measured quantities. As input for the simu- 
lation a spectral function is needed. Because the 
spectral function is best known at small E m and 
p m , one restricts this measurement to the IPSM 
region. As boundaries E m < 80 MeV and p m < 
300 MeV/c are chosen, the same as used in previ- 
ous analyses at SLAC and Jlab [12113114] . Then 
the nuclear transparency is obtained by the ratio 
of the measured and simulated yield integrated 
over the chosen (E m ,p m ) region: 

_ J v dp m dE m N exp (E m , p m ) 



T A (Q J 



J v d Pm dE m N** m (E m , Pm ) ' 



(4) 



In the simulation a spectral function has to be 
chosen which should correctly describes the ex- 
perimental distribution in the region of interest. 
In the previous analyses an IPSM spectral func- 
tion was used which factorizes into an E and k 
distribution. The ^-distribution is described by 
a Lorentzian and the momentum distribution is 
derived from a Wood-Saxon potential whose pa- 
rameters are adjusted to data measured at Saclay 
|17j . This spectral function does not account for 
any depletion of the shells due to N-N correla- 
tions. To take this into account the simulated 
yield for carbon is divided by a factor of e SRC 1.11 
± 0.03. This factor was used in all the previous 
analyses and was estimated from an early exami- 
nation of N-N correlations in 12 C and 16 |18I19| . 
The inverse of e c corresponds to the occupa- 
tion number, which modern many-body theories 
predict to be lower, w 80 - 85% |2I3) . This number 
is larger than the spectroscopic factor giving the 
occupation of a state; it contains the background 
strength, also caused by correlations, which can- 
not be attributed to a specific orbit |2()I21| . The 
above considerations emphasize that the number 
of nucleons missing due to absorption in the nu- 
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Figure 3. Momentum distribution in the region 
of 0.03 GeV < E rn < 0.08 GeV obtained from 
the data, the CBF theory (solid) and the IPSM 
(dashed). Three data points (circles) are from 
data focusing on the high p m region |5I22) (s. sec. 

19- 



cleus cannot be distinguished from the depletion 
of the single-particle orbits due to SRC. 

The factor e SRC introduced in the IPSM spec- 
tral function to account for the depletion due 
to N-N correlations is certainly an approxima- 
tion. Comparison to data for higher p m (but still 
smaller than 300 MeV/c) reveal that the IPSM 
spectral function is still missing strength at E m 
around 0.05 GeV. To improve upon this approach 
a spectral function containing SRC from the be- 
ginning is employed in Eq. Q . This spectral func- 
tion is composed of a part due to SRC which ac- 
counts for 22% of the total strength as calculated 
in Local Density Approximation (LDA), and the 
IPSM spectral function mentioned above but re- 
duced by a factor (1 - 0.22) to ensure normaliza- 
tion. This approach circumvents the application 
of the extra factor e SRC . In addition, much bet- 
ter agreement between measured and simulated 
E m spectra is achieved. Using this spectral func- 
tion the five data points for the nuclear trans- 
parency of 12 C (solid symbols) are obtained which 
are shown as a function of the proton kinetic 
energy in Fig. |3 together with previous results 
obtained at SLAC [T^l (circles) and Jlab |13I14| 



(squares and diamonds) . The error bars shown in 
the figure contain the statistical and systematic 
uncertainty but not the model-dependent error of 
4.7 %. This applies also to the data points of the 
previous works. Since the previous experiments 
were analyzed using the same assumption and in- 
gredients the model-dependent error is the same 
for them, while it is somewhat lower in the case 
of using the CBF spectral function. 

The solid lines drawn in Fig. [21 are the result 
of the correlated Glauber theory which takes 
Pauli Blocking, dispersion effects and N-N corre- 
lations into account. For comparison also results 
from previous experiments fH|T2fl3| for iron and 
gold are shown. For all three nuclei and large 
proton kinetic energy (> 1.5 GeV) the theory de- 
scribes the data well within the error bars. At low 
energy there is remarkable agreement between 
theory and the experimental results obtained us- 
ing the CBF spectral function. 

In fig. |21 the momentum distribution, mainly 
covering the S 1/2 state in 12 C (E m = 0.03 - 0.08 
GeV) is shown for the data (squares), the CBF 
theory and the IPSM. Obviously the IPSM mo- 
mentum distribution starts to fail above p m ~ 
200 MeV/c. It exceeds the data for small p m by 
about 20 % because no correction for SRC is ap- 
plied. The CBF theory combined with the spec- 
tral function from IPSM as described above is in 
good agreement with the data. No extra renor- 
malization factor is needed. 

4. RESULTS AT HIGH E m AND p m 

To obtain the spectral function from eq. 01 a 
transparency factor of 0.6 was used for 12 C. The 
results presented here are restricted to the data 
measured in parallel kinematics to keep the cor- 
rections on the distorted spectral function small. 
The region of the onset of the A resonance is 
clearly visible in the data and its contribution is 
reduced by a cut. The experimental spectral func- 
tion measured for 12 C is shown in fig. ^together 
with the theoretical result obtained in the SCGF 
approach for nuclear matter at a density compa- 
rable to 12 C J3]. At low momentum, around the 
Fermi momentum, good agreement is observed. 
Interestingly this is also the case for the spectral 
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Figure 4. Spectral function for the 12 C nucleus. 
Experimental result for several momenta above 
the Fermi momentum are shown as solid lines 
with error bars. The dashed-dotted lines repre- 
sent the SCGF nuclear matter spectral function 
at a density of p = 0.08 fm" 3 |5I25| . 



Figure 5. A comparison between the experimen- 
tal result at k = 410MeV/c (solid lines with error 
bars), the theoretical spectral functions by Bcn- 
har et. al. [22 (dashed line), the SCGF result 
(dashed-dotted line) and the second order Greens 
function approach 3 (dotted). 



functions of refs. |24I3| . At higher p m , however, 
there are significant differences as can be seen 
in fig. [5] The spectral function from CBF the- 
ory has its maximum at around p} n /(2M) which 
would be the case for two nucleons leaving the nu- 
cleus after a hard collision due to SRC. However, 
the data show more strength at lower E m and 
the maximum of the spectral function is shifted 
to smaller E m values. This might be due to LRC 
in finite nuclei which are not treated in the nu- 
clear matter calculation. Further tensor corre- 
lations which are known to be more important 
than central short range correlations might lead 
to a relocation of strength. Comparison of the 
two Green's function approaches shows a remark- 
able improvement of the self-consistent treatment 
compared to the one which contains only 2nd or- 
der diagrams. This indicate that higher order 
diagrams are important and more complicated 
configurations than a simple hard interaction be- 
tween two nucleons are involved. 

From the spectral function the momentum dis- 
tribution and the number of protons found in the 
(E m ,p m ) region covered by the experiment can 
be obtained by integration. The lower integra- 
tion limit is set to E m = 40 MeV to avoid contri- 



butions from the single-particle region and LRC 
contributions; the upper limit is adjusted in such 
a way to exclude the A resonance. The same 
cuts are employed for theory. In fig. the ex- 
perimental momentum distribution is compared 
to the three theories. The result from the SCGF 
approach agrees best with the data whereas the 
other two calculations slightly over- or under- 
shoot the data. However, one should not forget 
that due to the choice of the integration limits in 
E m the comparison depends somewhat on the dif- 
ferent shapes of the experimental and theoretical 
spectral functions. 

Integrating the spectral function over the 
(E m ,p m ) region covered by the experiment {p m 
— 0.23 - 0.67 GeV/c) and using the integration 
limits discussed above, one obtains the correlated 
strength. This strength represents the number of 
protons which are located outside the IPSM re- 
gion. The numbers obtained from the data and 
from three theories are quoted in tabled To ac- 
count for FSI a correction of -4 % has to applied 
to the experimental value in table HI f7!25j . Good 
agreement is found between data and the CBF 
theory as well as with the SCGF approach. 

Comparison of the results to the ones obtained 
for heavier nuclei (Al, Fe, Au) shows that the 
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Figure 6. Momentum distribution of the data 
(circles) compared to the theory of refs. [3] (dots), 
@] (solid) and |2U (dashed). The lower integra- 
tion limit is chosen as 40 MeV, the upper one to 
exclude the A resonance. 



Experiment 


0.61 ±0.06 


Greens function theory 


0.46 


CBF theory |2] 


0.64 


SCGF theory g| 


0.61 



Table 1 

Correlated strength (quoted in terms of the num- 
ber of protons in 12 C.) 



shape of the spectral function for C, Al, and Fe 
ist quite similar. For Au a larger contribution 
from the broader resonance region is obvious and 
the maximum of the spectral function is shifted 
to higher E m . The correlated strength for Al, Fe 
and Au is 1.05, 1.12 and 1.7 times the strength 
for C normalized to the same number of pro- 
tons. This increase cannot be solely explained 
by rescattering but MEC's have probably taken 
into account. Another contribution may be com- 
ing from the stronger tensor correlations in asym- 
metric nuclear matter |26I27| . 
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